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Abstract 

We examine the scattering of charged pions from the trinucleon system at a 
pion energy of 180 MeV. The motivation for this study is the structure seen in the 
experimental angular distribution of back-angle scattering for 7r + - 3 He and 7r _ - 3 H 
but for neither 7r~- 3 He nor 7r + - 3 H. We consider the addition of a double spin flip 
term to an optical model treatment and find that, though the contribution of this 
term is non-negligible at large angles for 7r + - 3 He and 7r - - 3 H, it does not reproduce 
the structure seen in the experiment. 
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1 Introduction 



The measurement of the differential cross section for the elastic scattering of tt from 
3 H and 3 He has recently been extended into the backward hemisphere by Matthews et 
al, [[|]. For all four cases, at an incident pion energy of 180 MeV, there is a slight dip 
in the cross section at ~ 130° lab (or 5.8 fm -2 momentum transfer), but for 7r~- 3 H and 
7r + - 3 He there is also a very distinctive rise in the cross section at larger angles. Moreover, 
the authors stated that, due to the resolution of the spectrometer, the dip may be even 
deeper and sharper. They compared experimental data to scattering models by Kamalov, 
Tiator, and Bennhold and by Gibbs and Gibson |J. For 7r~- 3 H and 7r + - 3 He both 
models agree reasonably well with the cross section up to 130°, but both fail to predict 
the distinctive rise in the cross section at large angles. For 7r + - 3 H and 7r~- 3 He both models 
give acceptable agreement at all angles. Figure [T] compares the experimental data [|T], |J 
(squares) to a theoretical curve from Ref. || (dashed line). 

One is thus led to consider mechanisms that would give a significant rise in the cross 
section in the backward hemisphere for both 7r~- 3 H and 7r + - 3 He, but not for 7r + - 3 H or 
7r~- 3 He. Similar rises in large-angle cross sections have been seen for 7r + scattering from 
12 C and 16 O In the three- nucleon case the dependence of this effect on the charge of 
the pion and on the target may give an indication of the cause. 

The pion kinetic energy, T„-=180 MeV, is near the peak of the P33 resonance for n- 
nucleon scattering. At this energy the elastic scattering cross sections are in the ratio 
a(7T + p)/a(7T + n) « a(TT~n) / 'a{-K'p) Pd 9. Therefore, we expect 7r-nucleus scattering at 
T„.=180 MeV to be dominated by the ir + -p and 7[~-n amplitudes. We use the convention 
that even-nucleon reactions refer to 7r + - 3 He and 7r~- 3 H reactions (which have two domi- 
nant scattering centers), and odd- nucleon reactions refer to 7r _ - 3 He and 7r + - 3 H reactions 
(which only have one dominant scattering center at this energy) . We also use the nomen- 
clature of Schiff where an "even" or an "odd" nucleon corresponds, respectively, to one 
of the like nucleons or to the unlike nucleon (the neutron in 3 He or the proton in 3 H). 

Qualitatively we can understand the first dip (around 90°) as a minimum in the spin- 
independent amplitude due to the p-wave dominance of the pion-nucleon phase shifts. For 
this intuitive view, consider the single-scattering impulse approximation, where the basic 
dependence of the amplitude is given by the sum of the relevant 7rN amplitudes multiplied 
by a form factor. For 7i~-nucleon scattering, the amplitude consists of two incoherent terms: 
a non-spin flip term, f(6), and a spin-flip term, g(0). Near the peak of the P33 resonance 
f(9) ~ cos(#) (plus an s-wave contribution which shifts the minimum slightly away from 
90°), and g{9) ~ sin(0). 

In the approximation that the trinucleon system exists in a pure s-state, (which we 
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use throughout this paper), the amplitude g{6) arises entirely from the interaction of the 
pion with the odd nucleon since, in first order, the spin-flip amplitudes cancel from the 
even nucleons. For 7r~- 3 He this interaction is strong, thus near 90° cm. the minima in 
\f{d)\ 2 is significantly filled in by \g(6)\ 2 . For the even-nucleon interaction 7r + - 3 He, n + -p 
has the largest amplitude and the interaction with the odd particle is weak so that the 
filling of the minimum is considerably less. An analogous argument follows for 7r ± - 3 H. 

It has been suggested [[JJ that the rise in the cross section at large angles might be due 
to the interaction of the incident pion with both of the like nucleons, flipping the spin 
of each so that the spin of the pair is conserved. That is, for 7r ± - 3 He, the incident pion 
sequentially induces a single spin flip of each proton, thereby leaving the final pair with 
spin 0. In the common optical model treatment, a potential is constructed from the single 
scattering amplitude of the type discussed previously, and then the nuclear amplitude is 
obtained from the solution with this potential in a wave equation. Thus a spin projection 
change is absent from this type of treatment and must be calculated separately. For this 
reason double spin flip (DSF) scattering is not directly included in current optical models. 
The DSF, whose amplitude is coherent with f(9), should clearly have more effect on the 
cross section of the even-nucleon reactions than of the odd-nucleon reactions. As each 
spin-flip amplitude preferentially leads to scattering around 90°, the two scatterings will 
lead to a forward-backward peaked angular distribution of the scattered particles. The 
forward part will likely be much smaller than the amplitude from the first order optical 
potential, but at large angles the two contributions may be comparable. 

Franco || has investigated multiple spin flip effects for 7r- 4 He scattering in the Glauber 
approximation. However, the Glauber approximation is not applicable in the large angle 
scattering region where the DSF is expected to be of significance. 

We treat the DSF as a second order correction to the scattering model of Ref. || 
and calculate it in the distorted wave impulse approximation (DWIA). This calculation is 
presented in Section |3], after briefly reviewing the scattering model of Ref. || in Section |2|. 
The results are discussed in Section f|, and conclusions are presented in Section |5[ 

2 Basic Scattering Model 

In the case of moderately heavy nuclei the optical model can be extended to include, by 
direct calculation, a number of effects including finite range and medium modification of 
the pion- nucleon scattering amplitudes ||. The medium corrections can be understood, 
in an approximate sense, in terms of a shift in the energy at which the phase shifts are to 
be evaluated and a transformation from the pion-nucleon to the pion-nucleus frame (the 
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so called "angle transform"). The effect of the finite range of the pion-nucleon interaction 
is included by assuming plausible forms for the off-shell t-matrix for the 7rN interaction. A 
number of studies of the pion-nucleon interaction have attempted to determine the ranges 
(there is a different range for each spin-isospin partial wave). For a recent analysis using 
local potentials see Ref. |T(J|]. These ranges can also be treated as phenomenological and 
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Another important element in pion scattering theory is due the possibility of true 
absorption of the pion, converting its mass into energy, as opposed to the usual "opti- 
cal model absorption" where the incident particle is simply removed from the beam by 
inelastic scattering. A number of attempts have been made to include this effect from 
fundamentals [|l|, [OJ. In the present calculations we adopt the method used in Ref. 
of including an imaginary term in the potential proportional to the square of the density. 
Approximate values of the parameter can be estimated from fits to heavier nuclei (see 
Refs. [0,0). 



The optical model contains a long known |15| correction due to the fact that, since 
the t-matrix is used to describe each individual scattering, no nucleon can be struck twice 
successively. Thus, in a multiple scattering picture, there can be A scatterings the first 
time but only A-l for each subsequent order. If all scatterings are of the same strength, 
this correction can be made by solving a wave equation with a potential having an overall 
strength of A-l, instead of A, and then multiplying the resulting t-matrix by the factor 
A/(A-1), (the KMT factor). For an optical model in which there are only three nucleons 
involved, this effect is of considerably greater importance than for a heavier nucleus where 
this factor is close to unity. We do not attempt to do anything beyond what was used in 
Ref. 0. The forward cross section is normally the most reliably calculated in scattering 
theories and it was noted in Ref. [pi that a treatment of this correction involving factors 
of the order described above was necessary to obtain agreement with the data at forward 
angles. 

We have carried out a search over variations of the scattering parameters used in 
Table I of Ref. ||. No substantial rise in the backward direction was seen. In order to 
focus on the DSF contribution to the cross section, we have fixed the values at: Energy 
shift=0 MeV, angle transform parameter=l, p 2 coefficient=8.8 fm 4 , and the off shell 
ranges (s and p) at 600 MeV/c. For a discussion of the corrections which come in for 
scattering of pions from few-nucleons systems see Ref. [18] . 

For the KMT factor, the calculated value (Eq. A4 of Ref. M) is used. We emphasize 
that a factor of this size is necessary in order to get agreement in the forward direction. If 
a factor of unity is used the cross section forward of 90° is shifted toward smaller angles. 
However, see the discussion in the conclusion on this point. 
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The spin-flip amplitude, g(6), is calculated in the distorted wave Born approximation, 
where the distortion is due to scattering from the even nucleons. The appropriate poten- 
tial to use in calculating these distorted waves is unclear. The KMT treatment for the 
elastic scattering from two particles indicates that the potential should contain a factor 
of 1/2. Using only the wave functions from this calculation, however, no correction to 
the amplitude is ever made. Thus it is not clear if a factor of 1/2, unity, or an interme- 
diate value should be used. The use of 1/2 gave favorable results when calculating the 
odd-nucleon interaction cross sections in Ref. [[J and we use the same value here. 

In calculating the DSF amplitude the distortion is taken to be from the odd nucleon 
and one of the even nucleons. We use a factor of unity here. This is consistent with the 
study of single charge exchange |16| (in which the distortion is also due to one strong and 
one weak interaction), for which the results were best when the full optical potential was 
used. 



3 Double Spin Flip Calculation 



Let k and k' be the pion's initial and final center of mass momenta, respectively, and 
|k| = |k'| = k, k • k' = cos#. Let the two even nucleons have coordinates r x and r 2 , 
spin operators cr 1 and cr 2 , and wave function x( r i, r 2)- The pion's initial and final 
distorted wave functions are then \I/( + )(k, ri) and \&(~)*(k', r 2 ). Assuming closure over 
the intermediate states and plane wave propagation of the pion in the intermediate state, 
the double scattering amplitude is given by 



F(k,k') 



2tt 2 



dr 1 dv 2 dqx*(ri,r a 



-jq-ri 



where the pion propagator is 



/ 2 (q,k 2 )*(-)*(k',r 2 ) 
/ 1 (k 1 ,q)vI/(+)(k,r 1 )] X (ri,r 2 ), 

P(q) ' 



3 iqr 2 



P{q) x 



q 2 — k 2 



the pion-nucleon spin flip amplitude is 

/i(qj.>Qa) = iAv(gi) v(q 2 ) (Ti ■ (qi x q 2 ) 
and the off-shell form factor is 



v(q) 



a 2 + k 2 
a 2 + q 2 



(1) 
(2) 

(3) 
(4) 
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The quantity a describes the range of the pion-nucleon interaction, which we take to be 
600 MeV/c, and A is derived from the 7rN phase shifts [[17]] • Since f\ is an operator in 
nucleon spin space, we must compute the expectation value of the operator 

O = (oi • k x X q) (cr 2 • q x k 2 ) . (5) 



For a pure singlet case [19| 

(O) s=0 = - (k x X q) • (q X k 2 ) = -q 2 [(k x ■ q) (k 2 • q) - k x • k 2 ] . (6) 

We write 

|x(ri,r 2 )| 2 = p{rx)p(r 2 ) (7) 
for the target in its ground state, so that the DSF scattering amplitude is now 

F(k,k') = ^1 / d ri dr 2 dciq 2 v 2 (q)P(q)e^ r2 p(r 2 )e~ 1 ^ p( ri ) x 

v{h) v(h) [(kx ■ q) (k 2 ■ q) - k x ■ k 2 ] ¥~>(k', r 2 ) < + >(k, n) . (8) 
As the momenta are operators on the distorted waves, we set 

k^-iVj, (9) 

and define 

V>(k, n) = v(h)¥ + \k, n) ; ^(k', r 2 ) = V (fc 2 )tf(->*(k\ r 2 ) . (10) 
We use a partial wave expansion for the pion wave function 

¥ + \k,n) =4 7 r^z i y im (k)y^(f 1 )^(A;r 1 ) , (11) 

so that, by use of a double Fourier transform and ort honor mality of the spherical har- 
monics, 

V>(k,r x ) = 1 f + *' a /d» e -**^gW(k,*) (12) 
(27r) a - Vi j 

= - (a + k J 4tt 2^ i Y lm (k) Y^fa) / dr dn ^— ,(13) 
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where we identify 

/, % 2 / 9 ,o\ f , , r 2 k 2 7/(«n) jj(Kr)ty](kr) , , 

^ = -( a 2 + k 2 ) drdn ^ 2 , 2 (14) 

drr 2 Gi(r,ri) ^i(kr) , (15) 
and where contour integration gives 

Gi(r,r) = -la +k ) dn — — (16) 

7T v ' jo or + /c 

= — a (a 2 + A; 2 ) /i| + ' ) (iar > ) jiiiotr^ (17) 
for r > the greater of r and r'. In a similar fashion 

^(k',r 2 ) =47rJ2(-i) l Y lm (k')YL(h)Mkr2) , (18) 

as #H*(k,r) = *W(-k,r). 

Thus Eq. |8| in the DWIA may be written 



F(k,k') = ^ / rfr 1 rfr 2 rfqg 2 t; 2 (g)P(g)e iq - r2 p(r 2 )e- 4qri p(r 1 ) x 

[Vx • V 2 - (q • Vx) (q • V a )] ^(k, r x ) ^(k', r 2 ) . (19) 

The DSF amplitude is twice the value given by Eq. |Tj|, as the process may proceed in 
two time orders. 

3.1 Evaluation: First Technique 

While Eq. [19| may be calculated directly, it is more easily evaluated by integration by 
parts on ri and r 2 . We find 

A 2 r 

F(k,k') = -— / rfr 1 c/r 2 rfqg 2 w 2 (g)P(g)^(k,r 1 )e~ iq - ri ^(k',r 2 )e iqr2 x 

p'(n) p'(r 2 ) [h ■ f 2 - (q • r x) (q ■ f 2 )] . (20) 
Let us consider each term in brackets separately, and label them Fx and F 2 respectively. 
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The first term is calculated by expanding the wave functions, exponentials, and f i • r 2 
in partial waves. We find 

Fi(K k') = ~ (4vr) 3 £ (2L + 1) (C^Z) 2 ^l(cos B) I^fa, r 2 ) , (21) 



Li 



where 



dr^rirlp^^Mkr^Mkr,)^^), (22) 



and 



(•OO 

T] luh (r 1 ,r 2 ) = J o dqq 4 v 2 (q)P(q)j h (qr 1 )j h (qr 2 ) 
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^e 3l< (krj h£\k r> ) + \a (o? + 3k 2 ) 3l< (zarj h£\ia r> ) 



+ (« 2 + ^) [r< f t< (iarj fy ( >W>) + j z< (iarj /i|+ ) '(iar > )] (23) 

provided < (Z> + 3) where r > is the greater of r± and r 2 and / > is r > 's respective index. 
To calculate the second term of Eq. EH, we write 



q ■ f = 4vr ^ Ji(?r) (C^?) 2 F iM (q) F L * M (r). 



(24) 



ILM 



Expanding the functions as before, we find 
A 2 



^(k, k') = — 2 (4vt) 3 £ (2L + 1) (C?£?J (C°°° 2 ) P L (cos9) I LM (n, r 2 ) . 



(25) 



Thus we have 



F(k, k') = 32ttA 2 £ (2L + 1) P L (cos 6) x 



000 



We may further reduce this expression with Racah algebra, 

F(k,k') = 327rA 2 E( 2 £+ 1)Pl(cos6) x 



(26) 
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{ -2(Q ° L ° 1 ) 2 (C 1 ° L ° 1 ) 2 W 1 , L+1 (r 1: 



+ 



+ 



(c? 




L L-l 



oo o 

L L+l 



^0 
°1 L L-l 



( r<o o o 

1^1 L L+l 



h,L,L~i,L-i(r 1 ,r 2 ) 



<L,L,L+1,L+1 



(n,r 2 ) } . 



(27) 



As the first two indices on the integral term, /, are the indices for the pion's wave function, 
the angular momentum of the pion is conserved as it must be for the elastic interaction 
with a spin-zero system. 



3.2 Evaluation: Second Technique 



Alternately, one may directly calculate Eq. |19[ This calculation is straightforward, but 
lengthy. We use the covariant components of the operator V in the spherical basis. For 

<7= (-1,0,1), 

V x ■ V 2 = £ (-IT V« V m a , (28) 



(q- Vi)(q- V a ) = 



r 000 r aa'0 v ( ~\ , 1 V~7 i \ M ^0 fiaa' M V ( *,\ 
WlO !! r 00iqj + —7= 2^\ l ) U 112°ll 2 r 2MiqJ 

v5 M 



^i-v a + J|E (-i) M c-'^y 2M (q)v« v« 



(29) 



Let i*s(k, k') be the contribution to -F(k, k') from Eq. ^8] and let Fo(k, k') be the contri- 
bution from the second term in Eq. |29|, so that 



F(k,k') = -F 5 (k,k') + F D (k,k'). 



(30) 



From Reference [EOfl we have 



1 A + l 
2A + 3 



A 



2A - 1 



CaTa-i ^aO" ») iA_i Ae+0 .(f i ) 



FiW s _ X_ Mkn) ; s A^n) + A±l ^ (31) 
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If we proceed in a manner similar to that in Section |3.1], we find 



and 



F s (k,k') = 327rAj^P i (cos^) f dr x dr 2 r\r\ p 2 (n) p 2 (r 2 ) x 
l J 

(' + !) Vi+u+Mi **) + 1 r 2 ) G?,(n) G,(rj,)] (32) 

dridr 2 rlrl p 2 ^^ p 2 (r 2 ) x 



F D (k,k') = -32vt^|aJ^(-1)'PKcos^) / 

[ (2/ + 3)(/ + l)77 i+im (r 1 ,r 2 )F,(ri)^(r 2 )Q° 1 ^ 1 2 { J T ^ } 
+/(2/-l)77 i _ M _ 1 (r 1 ,r 2 )G i (r 1 )GKr 2 )^° li 1 °{ J V ,1,} 

+ 2^/ (Z + 1) (2/ - 1) (21 + ^.^(n, r 2 ) G;( r i) F,(r a ) C^^° { « 'i 1 ( ^ }] (33) 
where 7y Zi Za (ri,r 2 ) is as previously defined. 

4 Results 

We performed calculations with the density of 3 He taken from the solution to the Faddeev 
Eqs. The two treatments presented in the previous section were evaluated and agree 
to within the expected numerical accuracy. The first method depends on the numerical 
calculation of the derivatives of the three-body density while the second does not. The 
derivative operators act on the pion wave function in the second method while the first 
method has a more direct form to calculate, i.e., it is expressed as an expectation value of 
a local operator. It is easier to see how to include the 5-function effect, discussed shortly, 
with the second method. 

Figure [I] shows the results of adding the DSF to the basic optical model for all four 
scattering cases. The dashed lines correspond to the optical model only and the solid lines 
include the DSF term. We see that the addition of the DSF term is indeed significant at 
large angles for the even nucleon cases. However, it does not give the structure seen in 
experiment. 

Figure |2] compares the scattering amplitude of the DSF to that of the basic optical 
model for 7r + - 3 He. The dashed line represents f(9) from the optical model only, the solid 
line the DSF amplitude, and the dash-dotted-line the DSF amplitude with the (^-function 
removed (discussed below) . We see that in the region of the first minimum both the real 
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and imaginary parts of the DSF amplitude are passing through zero, and thus have little 
effect in this region. 

One may speculate that a cancellation between the two mechanisms might produce a 
minimum at 130° and that the present calculation does not have the correct phase. To 
check this possibility, we introduced an arbitrary phase difference between the two terms, 
but still found no case which gave the characteristic minima. 

The asymmetries are shown in Fig. |3|. Since polarization phenomena are typically 
sensitive to small corrections, one might expect important corrections from the DSF; but 
we note that its addition has little effect since the asymmetry is only large around the 
minimum in the cross section where the DSF amplitude vanishes. 

It has long been known that there is a 5-function in the s-wave (nucleon-nucleon) 
portion of the one-pion-exchange potential (see e.g. Ref. |23[ for a discussion of this 
^-function and its removal). This 5-function also exists in the p-wave-p-wave part of pion 
double scattering ^|, especially visible in double charge exchange |TP|, |TJJ. 

If we refer back to Eq. |8], we see that the 5-function piece of the DSF amplitude results 



from the monopole portion of the expression [(ki • q) (k2 • q) — ki • k2] . In Section |T2 
we explicitly expanded this expression in terms of the V operator. Thus from the second 
treatment it is easy to see that the 5-function should be removed from the -Fs(k, k') term 
only. To do so, we make the replacement 



Q 2 



q 2 — k 2 



q 2 — k 2 



q 2 — k 2 



(34) 



which is equivalent to replacing q 4 with k 2 q 2 in the integral expression for Tj l j (ri, r%) in 
Eq. [23|. The term Fo(k, k') remains unchanged. 

As seen in Fig. ^|, this correction essentially reverses the sign of the real part of the 
scattering amplitude, while having a small effect on the imaginary part. However, it is the 
imaginary portion of the scattering amplitude that dominates the cross section at large 
angles where the DSF is of significance. Thus the ^-function correction has a minimal 
effect on the scattering cross section, decreasing it by less than 4 percent. In principle 
this correction should also be made in the basic optical model, akin to the Lorentz-Lorenz 



effect p5| , 27] at resonance energies. However, such a consideration is beyond the scope 
of this paper. 

It is known that the Faddeev densities, while representing the exact solution to the 
three-nucleon system expressed in terms of nucleon degrees of freedom, do not provide 
a completely accurate description of the electron-scattering cross sections. The problem 
is clearer at high momentum transfer near the first zero of the angular distribution |28| . 
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The difference is often ascribed to meson exchange currents [2^]. However, in order to 
explore the sensitivity of the results to the density used, we also performed calculations 
using the charge densities of Ref. JHJ corrected for the finite size of the proton using the 



proton charge parameterization from Ref. [31]. The two densities are compared in Fig. 



Only the proton density was changed, the neutron density remaining that of the Faddeev 
calculations. 

Figure [5] shows the comparison of the results of 7r + - 3 He scattering for the case in 
which the electron-scattering density was adjusted to have the same r.m.s. radius as 
the Faddeev density. The introduction of the electron scattering density at most shifts 
the cross section downward, and has little effect on the shape of the cross section at 
large scattering angles. For the calculation in which the two densities were used in their 
unmodified form the difference is even less. Thus it seems unlikely that the minimum is 
a form-factor effect, at least within the span of currently accepted functions. One can 
find a fit to the data by allowing an arbitrary form factor variation |32| but the density 
which results appears to be unphysical, having essentially no support for small values of 
the radial variable, i.e., it has a complete hole in the center. 



5 Conclusions 

We find that the double-spin-flip correction gives a non-negligible contribution to the large 
angle scattering cross section in the case of 7r + - 3 He. It should be included in any serious 
attempt to explain quantitatively the back-angle data on the three-nucleon system, and 
perhaps on other nuclei as well. However, because of the shape and phase of this coherent 
contribution, it does not appear that it can explain the dip and subsequent rise seen in 
the case of two strong scatters. 

As one possible direction for further work to explain the effect, we observe that if the 
KMT factor is set to unity, while the agreement of the forward cross section with data is 
completely ruined, there is a slight dip at the appropriate back angle. This suggests than 
an improvement to the basic optical model may lead to the resolution of this issue. An 



approach which may offer hope is that of Garcilazo | 33fl , but such a treatment involves a 
complete reformulation of the scattering theory. 

The single energy shift used in the present calculations might also cause problems. In 
Ref. || it was found that the shift in energy was dependent on the pion-nucleus angular 
momentum which is an alternate expression of the results of earlier work in which it was 
a function of momentum transfer [34|. 



These two observations lead one to think that the next appropriate step might be a 
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re-examination of the scattering theory to provide a more detailed representation of the 
physics. 
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especially the Theoretical Quarks, Hadrons, and Nuclei Group, for their hospitality. We 
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Figure 1: Comparison of experimental data from Ref. |] (squares) to a prediction from 
the optical model of Ref. P] with DSF term (solid line) and without DSF term (dashed 
line) . 
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Figure 2: Scattering amplitudes for 7r + - 3 He. The dashed line represents f(9) from the 
optical model, the solid line the double spin-flip mechanism (solid line), and the dash- 
dotted line the DSF with the 5-function removed. 
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Figure 3: Asymmetries for tt + 3 He and 7r~ 3 He scattering from the optical model only 
(dashed line) and with the DSF term (solid line). The data are from Refs. and J23L 
For 7r + there are additional data points at forward angles which are negative (not shown). 
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Figure 5: Comparison of the results of the optical model with the Faddeev density (solid 
line) and the electron scattering density (dashed line). The radius for the electron scat- 
tering density has been rescaled the same r.m.s. value as the Faddeev density. If no such 
rescaling is done, the difference is considerably less. 
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